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Abstract. The model of hydrogen migration and of the reactions of hydrogen atoms 
with electrically active impurity, developed earlier, has been applied to simulate hydrogen 
diffusion and passivation process during plasma deuteration of silicon substrates doped 
with boron. The calculated deuterium concentration profiles agree well in the length 
of the passivated region with the experimental data obtained on treatment in hydrogen 
plasma at a temperature of 200 °C for 5, 10, and 15 minutes. On the other hand, to 
achieve a good fit to the abruptness of the calculated profiles between the passivated 
and unpassivated regions, it is necessary to suppose that the values of the parameters 
that describe the absorption of hydrogen interstitials by electrically active dopant atoms 
decrease with increase in the depth of the passivated region. For example, nonuniform 
spatial distributions of nonequilibrium point defects generated during plasma treatment 
can lead to a spatial dependence of hydrogen absorption. 
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1 Introduction 

The increase in the cost of traditional energy resources has been a worldwide tendency dur- 
ing the past years. Therefore, investigation and implementation of alternative renewable 
energy sources are particularly urgent now. Solar cells provides such type of promising 
alternative energy sources. Here, new modules based on silicon layers constitute a re- 
liable, proven, sustainable, and environmentally friendly source of energy. It is worth 
noting that the best crystalline silicon photovoltaic modules are 5 % more efficient than 
the best modules based on polysilicon films [H |2] . It is especially important for solar cells 
used in outer space. The treatment of silicon photovoltaic layers in a hydrogen containing 
gas plasma results in a further increase in the solar energy-to-electricity conversion ratio. 
All of the above-mentioned factors demonstrate the importance of elucidating the role of 
hydrogen in the evolution of defect-impurity system in the near-surface region of silicon 
layers including passivation of dangling bonds, undesirable defects, and dopant atoms. 
The main goal of this work is to analyze theoretical models and carry out calculations 
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of hydrogen diffusion and hydrogen passivation of electrically active impurity atoms in 
silicon crystals doped with boron. 



2 Model 

The model of hydrogen diffusion with account for the passivation of electrically active 
impurity atoms was developed in [3l H] . It is worth noting that in contrast to other models 
of hydrogen diffusion (see, for example, [5]), in the model of Velichko et al. [HI [4], as well 
as in the model of Zhang [6], the values of hydrogen diffusivity are approximated from 
high-temperature experimental diffusivity data obtained in [7]. Thus, the contradiction 
between the diffusivity data obtained for low- and high-temperature treatments (see [6] 
and references therein) is eliminated. It is supposed in [3l H] that there are two fluxes of 
hydrogen atoms: "slow" diffusing hydrogen species, responsible for diffusion in the near- 
surface region with high a hydrogen concentration, and "fast" diffusing species, responsible 
for the hydrogen diffusion into the bulk of a semiconductor in the low-concentration region 
of the hydrogen concentration profile. It was supposed that the fast diffusion occurs due to 
the long-range migration of nonequilibrium interstitial hydrogen atoms in singly negatively 
(H~), neutral (H^), and singly positively (H"*") charge states. Due to the high mobility 
of electrons (holes), the mass action law is valid for conversions between different charge 
states of hydrogen interstitials. The long-range migration of hydrogen interstitials results 
in the supersaturation of the bulk of a semiconductor with hydrogen atoms. Therefore, 
it is supposed in [3l H] that these fast diffusing species are responsible for the passivation 
of electrically active boron atoms due to the following reactions: 

A- + H+^(AH)^ (1) 

A- + H^ ^ (AH)^ +e~, (2) 

where A~ is the acceptor atom in the substitutional position; (AH)^ is the neutral im- 
mobile "boron-hydrogen" complex, and e~ is the electron. Then, to calculate the total 
hydrogen concentration profile and concentration profile of the electrically active boron 
the following system of equations can be used: 

1) the conservation law for trapped hydrogen atoms: 

ac^^^ _ fc^^^(x) c^^^(x,t) 

+ kf^^k^^^^ix) C"^"" {x, t) + S^^{x, t) - G^^ix.t) , 

2) the conservation law for substitutionally dissolved dopant atoms: 

^^1^ = -kr^k^'^^\x)C C^^x (x, t) , (4) 

3) the stationary diffusion equation for nonequilibrium hydrogen interstitials: 
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d_ 
dx 



dx 



(5) 

- A;f^-^A:^^^^(x)CC^^^(x,t) + = 0, 

4 ) the diffusion equation for slow migrating hydrogen species: 

where 

^-1 I 1 I aHIS+ ^ 

\X) - pHis- + 1 + ' 

^"""(X) = ^^^^ . (9) 



Here C^"^ is the total concentration of hydrogen atoms; C^^^ is the total concentra- 
tion of the hydrogen atoms trapped by immobile sinks and boron atoms; C and Cb stand 
for the concentration of electrically active dopant atoms which undergo passivation and 
the concentration of charged species with the opposite type of conductivity, respectively; 
C^^ and C^^^ represent the total concentration of interstitial hydrogen atoms and the 
concentration of hydrogen interstitials in a neutral charge state, respectively; C^^ is the 
concentration of the slow diffusing hydrogen species; x is the concentration of charge 
carriers (holes p or electrons n for passivation of acceptor or donor impurities, respec- 
tively) normalized to the concentration of intrinsic charge carriers nf, S^^ and G^^ are 
respectively the rates of direct trapping and detrapping of hydrogen atoms introduced 
into the near surface region of silicon substrate by platelets or other extended defects; 
G^^ is the generation rate of noncquilibrium hydrogen interstitials in the surface region 
due to the plasma treatment and dissolution of platelets or other defects which incor- 
porate hydrogen atoms; G^^ is the generation rate of hydrogen atoms participating in 
slow diffusion; df^ and d^^'"{x) ^'^^ the diffusivity of hydrogen interstitials in an intrinsic 
semiconductor and normalized concentration dependence for diffusivity of this species in a 
doped semiconductor, respectively; D^^ is the diffusivity of the slow migrating hydrogen 
species; r/^^ = {kf^)~^ is the average lifetime of noncquilibrium hydrogen interstitials in 
an intrinsic semiconductor; kf^ and k^^^{x) ^"^^ the coefficient of absorption of hydrogen 
interstitials in an undoped semiconductor and the normalized concentration dependence of 
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this coefficient in a doped semiconductor, respectively; kf^^^ and k^^^^i^x) are the coeffi- 
cient of absorption of nonequihbrium hydrogen interstitials, when passivation of impurity 
atoms occurs in the near intrinsic sihcon, and the normahzed concentration dependence 
of this coefficient in a doped semiconductor, respectively. 

The empirical parameters and (3^^^ describe the relative contribution of neg- 

atively and positively charged hydrogen interstitials, respectively, to the total hydrogen 
diffusion in comparison with the contribution of neutral interstitials. The empirical pa- 
rameters /3^^^~ and (5^^^^ respectively describe the relative absorption of negatively and 
positively charged hydrogen interstitials due to unsaturated traps in comparison with the 
absorption of neutral interstitials. At last, the empirical parameter (3^^^ describes the 
relative absorption of positively charged hydrogen interstitials during the passivation of 
electrically active boron atoms in comparison with the absorption of neutral interstitials. 

It is important to note that the concentration dependences d^^'^{x)y k^^'^i.x)) and 
^HiAc^^ are smooth and monotone functions of x [8]. In addition, these functions have 
the form traditionally used for the presentation of effective diffusivity of substitutionally 
dissolved dopant atoms in processing simulation codes (see, for example, |i)J). Due to 
these features, the system of equations ([3]), (jl]), ([5]), and ([6]) is very convenient for nu- 
merical solution. It is also worth noting that the concentration of charge carriers x can 
be calculated either from the assumption of local charge neutrality (fTOj) or, more exactly, 
from the Poisson equation for electrostatic potential Lp. 



3 Simulation of hydrogen diffusion and boron passi- 
vation 

To illustrate the efficiency of the model developed in [3l H] for describing the hydrogen 
diffusion and passivation of boron atoms, the simulation results for experimental data of 
Tong et al. [10] are shown in Figs. [H El [3l HI and [51 In [10], deuteration was carried 
out in a theta-pinch plasma or a neutral atom gun. Deuterium profiling was carried out 
using SIMS with a cesium-ion source. The lower limit of the measurable concentration of 
deuterium was 10^ /xm~^. The deuteration temperature was chosen to be 200 °C and the 
surface treatment durations were 5, 10, and 15 minutes. 

To obtain the calculated deuteration concentration presented in Figs. [H El [HI HI and 
[5l we used an analytical solution of Eq. (El) for the case of the source provides permanent 
time-independent generation of a slow diffusing hydrogen species in the thin surface layer 
of silicon substrate. Thus, the flux of the slow component is approximated analytically. 
The reflecting boundary condition on the semiconductor surface is imposed on the low 
diffusing species with concentration C^^ . It can be seen from the experimental profiles 
that this type of a boundary condition is in agreement with the shape of hydrogen distri- 
bution in the near surface region. On the other hand, the finite-difference method [TT] is 
applied to find a numerical solution for the system of Eqs. ([3|), (jl]), and ([5]). It is supposed 
that generation of nonequihbrium interstitials occurs due to the plasma immersion ion 
implantation of hydrogen in the thin surface layer. In addition, the hydrogen interstitials 
can be generated due to the rearrangement or dissolution of platelets which incorporate 
hydrogen atoms. Unfortunately, the energy of hydrogen ions was not mentioned in flU] . 
However, the thickness of the layer where generation occurs is negligible in comparison 
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with the average migration length of hydrogen interstitials and character dimensions of 
the passivation region. Therefore, to obtain a numerical solution, the Dirichlet bound- 
ary conditions on the surface and in the bulk of the semiconductor are imposed on the 
rapidly diffusing neutral nonequilibrium deuterium interstitials with concentration C^^^ . 
Numerical computations were carried out on a ID simulation domain [0,xb], where xb 
and mesh point number were equal to 1.4 fim and 7001, respectively. Also, the time grid 
with 5000 equal steps was used to obtain a numerical solution. 

As can be seen from Fig. [H the calculated deuterium concentration profile is in good 
agreement with the experimental data of [10] . The following values of simulation param- 
eters were used to fit the calculated curve to the experimental deuterium concentration 
profile: Qf^ = 1.2x10^^ cm'^; Df^ = 2.0x10"^ /Jia^/s; df^ = 7.251 fim^/s; /f^ = 4.2 
/im; Cf^^ = 1.02x102 /im'^; = 0; /3^^+ = 1.0 xlO'^ Z?^^^" = 0; /3^^^+= 0; kf^^ 

= 0.21x10-3 /imVs; /^^^ = 0.4; C„= 1.3x10^ fim-^. Here Qf^ is the dose of hydro- 
gen atoms which are responsible for slow diffusing species; Cg^^ is the concentration 
of deuterium interstitials in the neutral charge state at the surface; If^ = ■>/ d^W^^^ is 
the average migration length of hydrogen interstitials in an intrinsic silicon; Cu is the 
concentration of uniformly distributed boron atoms in a substitutional position before 
hydrogenation. 




Figure 1: Calculated profiles of the total deuterium concentration (solid line) and unpas- 
sivated boron concentration (dotted line) for hydrogenation in the gas discharge plasma 
at a temperature of 200 °C for 5 min. The experimental data (black circles) are taken 
from Tong et al. [TO] . 

It follows from the fitting procedure that the increase of k^^"^ or (3^^^ results in 
a more abrupt deuterium profile between the passivated and unpassivated regions. The 
value of hydrogen diffusivity in intrinsic silicon df^ = 7.251 /im^/s for a temperature of 200 
°C was taken from Because the average migration length of deuterium interstitials 
in an intrinsic silicon l^^ is greater than 1 fim (it can be seen from Fig. [T] that the 
thickness of intrinsic silicon to the end of the passivation process is greater than 0.4 
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yum, and the thickness of hydrogenated layer is approximately equal to 0.6 /xm), we use 
in the simulation procedure the value If^ = 4.2 /xm. Then, the average lifetime of the 
nonequilibrium hydrogen interstitials in an intrinsic silicon r/^^ is equal to 2.43 s, i.e., 
r/^^ is significantly smaller than the duration of hydrogenation (5 min). This value of 
r/^^ confirms the correctness of the assumption about the quasistationary distribution of 
hydrogen interstitials due to their high mobility. 

In Fig. [21 the calculated deuterium concentration profile after thermal treatment for 
15 minutes is presented. The same values for diffusion of hydrogen interstitials and 
the same passivation parameters were used in these calculations. As can be seen from 
Fig. [21 the calculated deuterium concentration profile is in good agreement with the 
experimental data of [10], as concerns the length of the passivated region. On the other 
hand, the abruptness of the calculated profile at the boundary between the passivated 
and unpassivated regions is very high and disagrees with the experimental data. Similar 
calculations of the deuteration distribution for 10 minutes show that in this case the 
profile abruptness is also greater in comparison with the experimental one. 




Figure 2: Calculated profiles of the total deuterium concentration (solid line) and unpas- 
sivated boron concentration (dotted line) for hydrogenation in the gas discharge plasma 
at a temperature of 200 °C for 15 min. The experimental data (black circles) are taken 
from Tong et al. [TO] . 

It is rather difficult to explain the obtained disagreement at the boundary between 
the passivated and unpassivated regions after plasma treatment for 10 and 15 minutes. 
Indeed, to reduce the abruptness of the calculated deuterium profile, one can decrease the 
value of kf^^ and/or (3^^^. In Figs. [21 [H and [5l the calculated deuterium concentration 
profiles with k^^"^ = 0.2x10"^ fim^/s and (3^^^ = 0.12 a. u. after thermal treatment for 
5, 10, and 15 minutes are presented. 

It can be seen from Figs. [3l [H and [5] that the deuterium concentration profile 
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Figure 3: Calculated profiles of the total deuterium concentration (solid line) and unpas- 
sivated boron concentration (dotted line) for hydrogenation in the gas discharge plasma 
at a temperature of 200 °C for 5 min. The decreased values of kj^^^ = 0.2x10"^ /um^/s 
and = 0.12 a. u. are used in simulation of hydrogen diffusion. The experimental 

data (black circles) are taken from Tong et al. [TO] . 



calculated for the decreased values of kf and (3 agree well with the experimental 
data if plasma treatment lasts for 10 and 15 minutes whereas for the 5-minutes duration 
the slope of the experimental profile is steeper. At present, we do not know the physical 
reasons for the decrease in kf^^ and (3^^^ with increase in the duration of deuteration 
and respectively in the depth of the passivated region. Indeed, the reactions between the 
hydrogen and boron atoms are independent of the distance from the surface. Perhaps, 
there is an error in the measurements of the deuterium concentration profile, or the 
nonequilibrium point defects generated in the near surface region due to plasma treatment 
change the conditions for passivation. For example, generated silicon self-interstitials can 
kick-out boron atoms from their substitutional position (Watkins effect). This problem 
requires a further investigation. 

It is worth noting that simulation of the experimental data of Tong et al. pXlJ was 
also carried out in |6j. The model of hydrogen diffusion proposed in [6] also takes into 
account two different fluxes of hydrogen species. It is supposed that "fast" diffusion occurs 
due to the migration of atomic hydrogen in different charge states H~, H^, and H"*". To 
describe "slow" diffusion in the high concentration region, the mechanism of formation and 
migration of hydrogen complexes with mobile traps generated during plasma treatment 
is used. It is supposed in [6] that silicon vacancies trap hydrogen atoms and form mobile 
complexes V-H, which are responsible for the "slow" diffusion. In contrast to the present 
paper, only one quasichemical reaction, i.e.. 



7 




Figure 4: Calculated profiles of the total deuterium concentration (solid line) and unpas- 
sivated boron concentration (dotted line) for hydrogenation in the gas discharge plasma 
at a temperature of 200 °C for 10 min. The decreased values of kf^^^ = 0.2x10"^ ^m^/s 
and f]^^^ = 0.12 a. u. are used in simulation of hydrogen diffusion. The experimental 
data (black circles) are taken from Tong et al. [TU] . 



A- + H+^(AH)\ (12) 

is used to describe passivation of boron atoms. It is also supposed in [6] that reaction 
f fT2|) as distinct from reaction ([1]) is reversible, i. e., the detrapping of hydrogen atoms 
can occur. On the other hand, it was found in [12] that the dissociation energy of the 
(BH) complex is equal to 1.28 ±0.03 eV, i.e., it is high enough in comparison with the 
value of ksT = 0.04077 eV for a temperature of 200 °C. Here ks is the Boltzmann 
constant. Therefore, in reactions ([1]) and ([2]) the process of detrapping is omitted. As can 
be seen from the paper of Zhang [6], a full fitting of calculated deuterium concentration 
profiles to the experimental ones is not achieved, because the abruptness of the calculated 
deuterium distributions at the boundary between the passivated and unpassivated regions 
also disagrees with the experimental data. We hope that new, more precise measurements 
of deuterium distributions in silicon substrates doped with boron can solve this problem. 



4 Conclusions 

The model of hydrogen migration and reactions of hydrogen atoms with electrically active 
impurity which had been developed in [3l H] was applied for simulating the hydrogen 
diffusion and passivation process during plasma deuteration of silicon substrates uniformly 
doped with boron. For comparison, the experimental data of Tong et al. for plasma 
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Figure 5: Calculated profiles of the total deuterium concentration (solid line) and unpas- 
sivated boron concentration (dotted line) for hydrogenation in the gas discharge plasma 
at a temperature of 200 °C for 15 min. The decreased values of kf^^^ = 0.2x10"^ /im^/s 
and 13^^^ = 0.12 a. u. are used in simulation of hydrogen diffusion. The experimental 
data (black circles) are taken from Tong et al. jlOj . 

treatment at a temperature of 200 °C for 5, 10, and 15 minutes was used. The calculated 
deuterium concentration profiles agree well with the experimental data for the entire 
time of plasma treatment if we suppose that the coefficient of absorption of hydrogen 
interstitials in the near intrinsic semiconductor kf^^^ as well the empirical parameter /3^^^, 
which describes the relative absorption of positively charged hydrogen interstitials in 
comparison with the absorption of neutral interstitials during the passivation of electrically 
active dopant atoms, are decreased with increase in the duration of deuteration and 
accordingly in the depth of the passivated region. It is possible that this decreasing 
occurs due to the nonuniform spatial distribution of nonequilibrium point defects, which 
are generated during plasma treatment and can influence the passivation process. If the 
invariable values of the simulation parameters are used for different durations of plasma 
treatment, the calculated deuterium concentration profile is in good agreement with the 
experimental data of (TU], as concerns the length of the passivated region. However, the 
abruptness of the deuterium concentration profiles at the boundary between the passivated 
and unpassivated regions disagrees with the experimental data either with treatment for 
5 minutes or 10 and 15 minutes. 
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